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We show that string theory with Dirichlet boundaries is equivalent to string theory 
containing surfaces with certain singular points. Surface curvature is singular at these 
points. A singular point is resolved in conformal coordinates to a circle with Dirichlet 
boundary conditions. We also show that moduli parameters of singular surfaces coincide 
with those of smooth surfaces with boundaries. Singular surfaces with saddle points indeed 
arise in the strong coupling expansion in lattice QCD. The kind of saddle point, which 
may be the origin of a singular point we need, is of infinite order. 
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Green proposed to study string theory propagating with surfaces bordered by Dirichlet 
boundaries, in order to reproduce the high energy behavior of scattering amplitudes in QCD 
0. Earlier references in the subject include those in ||. For a recent approach to this 
string theory, see @ . Is there any evidence based upon more theoretical grounds in support 
of the Dirichlet string theory, other than merely requiring account for "phenomenological" 
high energy results? The purpose of this note is to provide such evidence. 

Starting with QCD, either from 't Hooft's 1/N expansion [|J or from Wilson's strong 
coupling expansion on the lattice 0, one arrives at the conclusion that QCD may be 
described by a non-interacting string theory at the large N limit. Unfortunately this 
turned out not quite to be the case 0, at least some modifications are needed. Kazakov 
proposed a modified strong coupling expansion in [[7], in which not only more degrees of 
freedom are required, but also some singular saddle points appear on random surfaces. 
This approach was further developed in ||, for a recent review, see @. Despite subtleties 
in going from the strong coupling regime to the weak coupling regime, those saddle points 
nevertheless survive. It is these singular points which interest us here. We shall argue that 
these points could serve as the origin of Dirichlet boundaries. 

In the strong coupling expansion, a saddle point arises when 2n (n > 1) plaquettes 
are contracted cyclically on a link. The deficit angle at this point is 2n(n — 1). If this 
point survives the continuum limit, this deficit angle will induce a singular curvature, i.e., 
the curvature blows up approaching this point. Associated to a saddle point of a given 
order, there is a coefficient depending on the type of the saddle point. For the simplest 
saddle point with n = 2, there is an additional weight factor 1/N. There are higher order 
corrections in terms of 1/N. These can be explained as the coalescing of several saddle 
points of the same order. Now let us assume that this singularity can be described by a 
singular intrinsic curvature on the continuum surface. For a finite n, the curvature behaves 
as a 6- function centering at the saddle point. This kind of singularity will turn out to be 
too mild to induce a Dirichlet boundary. Nevertheless let us discuss this singularity first. 
Note that a saddle point of finite order is precisely a branch point of the same order found 
in the QCD2 string theory |TD]. To see this, consider a metric ds 2 = \z\ 2 ( n ~ 1 * ) \dz\ 2 . The 
curvature of this metric behaves as y/gR ~ S 2 (z). Use coordinates w = z n instead, the 
metric becomes smooth ds 2 = \dw\ 2 . The singular property is encoded in the fact that the 
z plane covers the w plane n times, and z = is a branch point of order n. It is instructive 
to define the deficit angle by 
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where C is the contour with a constant radius r. The first term is just the ratio of 
the circumference to the radius. In polar coordinates z = re 1 ^ , the metric reads ds 2 = 

is easy to check that our definition of the deficit angle yields a value 
27r(n — 1) in this case. To study an appropriate conformal field theory, which is needed 
in a string theory, on a closed surface with such a branch point, one may either resolve 
this mild singularity by going to a higher genus (as done in the QCD2 string theory), or 



one may simply insert an appropriate vertex operator [11]] at this point and deal with a 
smooth surface (in coordinates w). 

Now the kind of singular point inducing a Dirichlet boundary is different from a saddle 
point of any finite order. The reason is that the curvature blows up badly at the singular 
point. Without loss of generality, let us consider a metric 

ds 2 = r 2a d(j) 2 + dr 2 . (2) 

This metric is singular at r = when a < 1. We also require a > 0, otherwise r = 
is not really a point. A simple calculation yields R = 2a(l — a)r~ 2 . This curvature 
is so singular that the integral J ^fgK diverges at r = 0. Note that we have to put the 
restriction a < 1, otherwise J y/gR converges. Thus, this singular point can not be resolved 
by going to higher genus. This is not surprising, since the deficit angle as defined in (TJ) 
is 27r(r a_1 — 1), which approaches infinity as r approaches zero. If we are to explain this 
singularity as originating from saddle points, it must consist of infinitely many saddle 
points on the same link with arbitrary large order. We conjecture that this is the case, 
although it is hard to have a derivation directly from the strong coupling expansion. After 
all the strong coupling expansion has not been proven to have a continuum limit yet. 

We shall not assume any particular conformal field theory (CFT) in the QCD string 
theory here. To study a CFT with geometry (|!|), it is necessary to find a coordinate system 
in which the metric reads exp((p)\dz\ 2 . Then we need to study this CFT in the geometry 
\dz\ 2 by a Weyl rescaling. It is straightforward to see that the metric (||) can be written 
as r 2a \dw\ 2 , where w = y^r 1- " + i<p. Now w furnishes complex coordinates on the semi- 
infinite cylinder, and the boundary of this cylinder is at r = 0. Notice that if a = 1, 
one can not use these coordinates. The semi-infinite cylinder can be further mapped to 
the complex plane with a unit disk removed, by z = exp(w). The metric reads in these 
coordinates 



The Weyl factor exp(0) vanishes on the boundary \z\ = 1. The Liouville action, resulting 
from a Weyl rescaling, is infinite with this Weyl factor. We have seen that in a conformal 
coordinate system, the singular point necessarily blows up to a circle. If the CFT to be 
studied includes spacetime coordinates as a subset of conformal fields, the boundary 
conditions for these fields are Dirichlet, as this boundary represents just a previous singular 
point of curvature. We also note that the resulting conformal geometry is independent of 
parameter a. Formally, the Liouville action depends on a, however, we shall assume that 
this Liouville action be eventually canceled out by that from contribution of b-c ghosts. 
We postpone discussion of boundary conditions for these ghosts to a little later. 

The above consideration of resolving a singular point can be generalized to the sphere 
with a singular point, and in general to a sphere with many such singular points. An 
interesting question then arises. If any such singular geometry is resolved resulting in a 
sphere with boundaries, are there enough configurations of singular geometries to cover 
the moduli space of the conformal sphere with boundaries? The answer to this question 
is positive. We consider the sphere with two singular points as an example. Assume that 
r = and r = oo be two singular points on the sphere, in polar coordinates. The metric 
is approximately as in (|2|) around r = 0, with parameter a\. The metric is approximately 
r' 2a2 d(p 2 + dr' 2 around r = 1/r' = oo. a\ are all positive and less than one. We can write 
down such a metric readily 

ds 2 = (1 + r 2 )" 2 ((r 2ai + r 4 " 2a2 )# 2 + dr 2 ) , (3) 

r 2ai dominates r 4_2a2 when r — > 0; while r 4_2a2 dominates r 2ai when r — ► oo. The above 
metric is written in conformal coordinates as 



\dz 
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ds 1 = (l + r 2 )" 2 (r 2ai +r 4 " 2a2 )^-, (4) 



where 

r r dr 
z = e R+i( t>, R(r): 1 



( r 2ai _|_ r 4-2a 2 )l/2 ' 

Apparently R is a well-defined function of r. The singular point r = is mapped to \z\ = 1, 
the unit circle. The singular point r = oo is mapped to the circle \z\ = exp(R ) > 1. R 
is finite, and given by 

_ f°° dr 1 l- fll l-a 2 \ 



Q ( r 2ax +r 4-2a 2 )i/2 0F(4 - 2a x - 2a 2 ) v 4-2ai-2a2 4 - 2a x - 2a 2 

(5) 
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We thus see that the sphere with two singular points is mapped to a conformal sphere 
with two disks removed, and it is just the annulus. The single real moduli parameter is the 
radius of the outer boundary exp(i? ) 5 which is finite. The moduli parameter depends on 
both a\ and a^. Adjusting Oj, Rq can be sent to infinity. However, R can never be zero, 
so geometries as described by (^]) can never cover the whole moduli space. To cover the 
whole moduli space, one can replace the factor r 2ai -\-r 4 ~ 2a2 in (|3|) by (r 2ai +r 4_2a2 )p 2 (r). 
Here p(r) is a function with p(0) = p(oo) = 1. Now Rq is not determined by ([5]) but by 

dr 

It is possible to adjust p(r) to let Rq approach zero. This is done by increasing p(r) for 
all r except r = 0, oo. Geometrically, the result of doing this is that the circumference of 
a circle with constant r is increased. Notice that the distance between two points r = 
and r = oo remains the same. This picture then explains why in conformal geometry, the 
outer boundary becomes closer and closer to the inner boundary. The whole moduli space 
is covered in this way. The moduli parameter then depends not only on parameters dj, but 
also on details of geometry. However, the existence of boundaries is due to the existence 
of singular points. 

These considerations generalize to a sphere with many singular points. It is hard and 
unnecessary to work out details of mapping singular geometries to conformal geometries 
with boundaries. We shall only argue that the whole moduli space of the latter is covered by 
varying singular geometries. Suppose there are N > 3 singular points, therefore there are 
N parameters ai characterizing these singularities. Together with the positions of these 
singularities, there are 3iV real parameters available. Since there is a global conformal 
invariance group SL(2,C) on the sphere, only 2N — 6 parameters out of 2N position 
parameters are relevant. So we are left with 3iV — 6 real parameters. This number is 
precisely the dimension of the moduli space of the sphere with iV components of boundary. 
As in the N = 2 case, moduli parameters are complicated functions of the parameters 
of singular geometries, and it may be impossible to cover the whole moduli space if we 
consider only one class of geometries, such as that in ([|). The moduli space can be covered 
by considering all possible singular geometries, again as in the N = 2 case. We believe 
that the conclusion presented in this section applies equally to Riemann surfaces of higher 
genus with singular points. 

After having shown that Dirichlet string theory is equivalent to a closed string theory 
with singular surfaces, we come to the issue of b-c ghosts. If the boundary conditions for 
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scalar fields are Dirichlet, what are the boundary conditions for b-cl Apparently, we can 
not impose Dirichlet boundary conditions on b-c consistently. These fields are fermionic, 
so Dirichlet boundary conditions can either be b = or c = on the boundary. If b = 0, 
then there is no nontrivial solution to the equation of motion db = 0, this would contradict 
the existence of moduli parameters. If c = 0, the trivial solution to dc = contradicts the 
existence of nontrivial conformal Killing vectors on, say, the disk. So we are left only with 
the possibility of Neumann type boundary conditions, namely cd = cd and b(dz) 2 = b(dz) 2 
on the boundary. We now show how this is consistent with consideration of singular points. 
Consider the ghost fields in the geometry (||). For the c field, there are two components 
(c r , d 3 ) in polar coordinates. In conformal coordinates (z, z), the two components become 
c z = zr~ a c r + izd and c z = zr~ a c r — izd . So both c 2 and c 2 become singular on the 
boundary, if c r ^ 0. If we demand that c r approach zero faster than r a (which approaches 
zero more slowly than r), we find the desirable Neumann boundary conditions. As for the 
b field, we find that two terms b rr (dr /dz) 2 and b r( j,(dr /dz)(d<p/dz) in b zz are vanishing 
when r — > 0, for dr/dz approaches zero. Neumann boundary conditions for this field then 
follow. 

We showed in that the Liouville action resulting from a Dirichlet boundary is 
essentially the same as the one from a Neumann boundary. In other words, Dirichlet 
boundary conditions do not change the central charge of the CFT. With a b-c ghost system 
obeying Neumann boundary conditions, it is possible to cancel the Liouville action from 
the matter CFT part. Together with our discussion above on the singular surfaces, we 
conclude that study of a string theory with such surfaces is made possible by studying 
the Dirichlet string theory. Because the latter makes perfect sense, it is then a reasonable 
hypothesis that singular surfaces contribute significantly in certain circumstances, such as 
in the QCD string theory. 

As argued by Green, as well as in || for a general case, the Dirichlet boundaries 
provide us with a mechanism to generate the appropriate high energy behavior observed 
in strong interactions, which an ordinary string theory fails to do. Why is this so from 
the standpoint of the strong coupling expansion? The rest of this note is devoted to a 
plausible answer to this question. 

First, let us remind ourselves how scattering amplitudes at high energies can be calcu- 
lated in the pure glue theory. Color singlet states are supposed to be glueballs. A glueball 
state, as an eigen-state of four momentum, is a superposition of many states with differ- 
ent numbers of gluons. Determination of the wave function of a particular glueball state 
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is a nonperturbative problem. Fortunately, as a high energy process at a wide angle is 
concerned, there must be a factorization theorem, similar to the one for meson scattering 
TIJ. This theorem states that the scattering amplitude factorizes into two factors. One 



of them has to do with the structure of the glueball, which may be called the distribution 
amplitude; another is the scattering amplitude between two individual gluons. The latter 
can be calculated perturbatively. Even the dependence of the distribution amplitude on 
the momentum transfer can be calculated perturbatively. The dependence of the whole 
amplitude on the large momentum transfer is largely determined by the two gluon scatter- 
ing amplitude. This is translated into the short distance behavior of the two point function 
of the gauge potential (A (l {x\)A v [x2)), which is calculated perturbatively. 

Instead of considering high energy scattering amplitudes of glueballs, one can consider 
correlation functions of gauge invariant local operators constructed from field strength. 
High energy behavior is translated into certain short distance behavior of these correlation 
functions. As far as any short distance behavior is concerned, one can use the operator 
product expansion, again calculated perturbatively []T3 



In the modified strong coupling expansion of Kazakov 0, a Lagrangian multiplier 
a, a hermitian matrix, is introduced to enforce unitarity of the link variable. Weight of 
each saddle point is the expectation value of a negative momentum of a, namely W n = 
(-^tra -2n ), where n is the order of the saddle point. It can be shown in the one plaquette 
model that solutions of the weight exist in both the strong coupling phase and the weak 
coupling phase . We are certainly interested in the weak coupling phase. Here one finds 
that for large /3, the inverse of the coupling constant, W n goes to (3~ 2n . This factor cancels 
the factor (3 2n from area of In plaquettes contracted on the link. This shows that all saddle 
points of higher order contribute in the weak coupling phase. 

Similar arguments lead to demonstration of the existence of the weak coupling phase 
in QCD4 ||, the physical model we are mostly interested in. Curiously, Kazakov showed 
that the weight W n is associated to local condensation of the gluon field. Absorbing a 
factor P 2n into the weight, he finds 

Wn = Pi^trFp + . . . (6) 

The quantity on the r.h.s. of (H) can be calculated perturbatively in the weak coupling 
regime, using the standard Feynman rules. Now consider correlation functions of gauge 
invariant local operators which in turn are constructed from plaquette variables, in the 
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framework of strong coupling expansion. Saddle points appear in any such correlation 
function, and the weight in the weak coupling phase is given by eq.(^). When two operators 
in the correlation approach each other, we expect that saddle points near these operators 
play a more and more important role. It is then expected that the operator product 
expansion in the lattice context is dominated by saddle points, this is consistent with the 
fact that the weight in eq.(^) is determined perturbatively. 

A lesson about the underlying CFT in the QCD string theory is readily drawn from the 
above discussion. In addition to the power-law behavior in the coefficients of the operator 
product expansion, there are logarithmic corrections, as functions of ln(p 2 / Aq CD ). The 
QCD scale Aqcd is the only free parameter in QCD. This parameter is replaced by the 
string tension T ~ m the string theory. So it is expected that similar factors 

involving ln(p 2 /T) appear in the string calculations, even when only one singular point on 
the surface is present. With the choice of the simplest CFT described by free fields X^ 1 on 
the world sheet, it is shown in [1]] and |J that no such logarithmic factors are obtained. 
Thus, some modified CFT is required to yield logarithmic corrections. This is in contrast 
to the belief that logarithmic corrections arise from renormalizing infinities in the string 
theory. If this were true, one would be forced to introduce another independent scale. It 
is questionable that such a scale can be related to the string tension in any reasonable 
manner. From the strong coupling expansion, it is obvious that the free CFT is not the 
right theory, for example back-tracking configurations are not suppressed in this theory. 
It is likely that the correct CFT for the QCD string, if any, must involve interacting fields 
A M and additional fermionic fields to give logarithmic corrections as well as to suppress 
back-tracking. 

In conclusion, no matter what the CFT may look like in the QCD string theory, 
singular surfaces must be taken into account in order to produce power-law behavior in 
high energy scattering amplitudes. We have shown that singular points can be resolved by 
introducing Dirichlet boundaries, so the theory can be studied systematically. 
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